Cosmological constraints for an Eddington-Born-Infeld field 
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We consider the Eddington-Born-Infeld (EBI) model here without assuming any cosmological 
constant. The EBI scalar field is supposed to play a role of both dark matter and dark energy. 
Different eras in cosmology are reconstructed for the model. A comparison is drawn with ACDM 
model using Supernova la, WMAP7 and BAO data. It seems that the EBI field in this form does 
not give good fit to observational data in comparison to the ACDM model. 



I. INTRODUCTION 



The cosmic acceleration is now considered to be one of the frontier quest of fundamental physics. Confirmed 
by observations ([1] and [2]), understanding of acceleration is yet to be satisfied in the regime of standard general 
relativity. Attempts to explain the acceleration could be performed by adding extra components of fluid into energy- 
momentum tensor part of the Einstein field equation. This extra component, dubbed dark energy, gives negative 
pressure so that it is able to drive the acceleration, see e.g. references in [3] for reviews. At smaller scales, a problem 
of an extra attractive gravity in galaxies and galaxy clusters shows up. Effects of extra gravity such as flat galactic 
• rotational curve, gravitational lensing, bulk velocity and structure formation are explained with dark matter [4]. On 
O" 1 observation side, the simplest model of dark matter and dark energy-the ACDM model is favored however suffering 
from fine-tuning problem. At present, nature of dark sectors is still unknown. There is another more radical way of 
acquiring acceleration. That is to modify the gravity term in the action (see [5] for recent reviews) . 

It is possible that dark energy and dark matter is only a single entity and it could effectively have different behaviors 
' at early and late times. This unified scenario is considered in Chaplygin gas model [6]. The Chaplygin gas model 
, however can not satisfy cosmic microwave background power spectrum and structure formation [7]. Another idea of 
unifying dark sectors was proposed recently by Banados [8] and [9]. The model called Eddington-Born-Infeld (EBI) 
gravity can account for both dark matter and dark energy components without additional degrees of freedom in 
energy-momentum tensor. In the model, Einstein gravity couples to Born-Infeld theory giving rise to a bi-metric 
theory. The second metric of the theory is generated from the Born-Infeld Christoffel symbol, which is solely 
responsible for dark sectors. The theory predicts a dust-like effective equation of state at large scale, while at late 
time it behaves like a cosmological constant. The theory can also accommodate fiat galactic disk-rotational curves. 
The model is motivated by a combination of Eddington's idea of purely affine theory of gravity without using metric 
[10], Born-Infeld- Einstein action [11], and the idea of magnetic spin symmetry breaking in presence of an external 
magnetic field. Considering that topological manifold is invariant under full diffeomorphism group of transformation, 
Ricmannian manifolds are invariant under smaller class of subgroup of metric isometries. The Eddington action which 
is diffeomorphism invariant hence is considered as an unbroken state theory [12]. Moreover, it is a ghost-free theory. 
Introducing <? M „ ^ to the gravitational theory would break this symmetry similar to having external magnetic field 
applying to a random spin system. The external magnetic field also breaks the symmetry of the spin system. If we 
let the metric couple to the Eddington action, the result is the EBI action. In the action, there is the Einstein-Hilbert 
part and the EBI part (see [8, 9, 13] for more detailed discussion). In context of anisotropic universe with Bianchi type 
I model, at late time, the EBI gravity effectively behaves like Einstein-Hilbert cosmology plus cosmological constant. 
The EBI term is stable at dark matter phase but also gives rise to anisotropic pressure and the perturbation decays 
oscillatory in time which differs from standard exponential decay case [14]. Considering the dark energy phase, the 
Born-Infeld as dark energy is not stable and it produces a very strong Integrated Sachs- Wolfe effect on large scales. 
This suggests that for the model to be viable, cosmological constant is needed in the action [15]. When adding a 
cosmological constant into the model, the model still predicts too large CMB fluctuation compared to WMAP5 data. 
However while restricting the EBI field as dark matter, the EBI model is a best fit with the ACDM prediction [16]. 
The idea that Eddington action is a starting point for general relativity is pursued further more when considering a 
Born-Infeld part and a cosmological constant but without having Einstein-Hilbert term. The idea was investigated 
in Palatini formulation to include matter fields. For homogeneous and isotropic space-time, the model present a 
non-singular cosmology at early time as well as non-singular collapsing of compact objects [17-20]. In such a scenario, 
Poisson equation is modified and Jean length is equal to the fundamental length of the theory. Also the critical mass 
for a blackhole to form is equal to the fundamental mass of the theory [21, 22]. 
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Here in this paper, we consider EBI model without cosmological constant as originally proposed in [9]. In fact, we 
believe that introducing a cosmological constant would make this model less attractive, as the model was introduced as 
a way to explain both dark energy and dark matter at the same time. In other words, adding a cosmological constant 
by hand would mean that the model achieves only half of the original goals it was introduced for. We re-emphasize 
the inviability of the original EBI model by fitting it with WMAP7, BAO and Supernova type la data. Compared to 
the study of [16], where the authors studied the growth of structure for these models (they studied the evolution of 
the cosmological perturbations during radiation and matter domination) we perform a study of the background and 
look for the constraints on it coming from the most recent data. In section II, we briefly describe the EBI model as 
a bi-metric theory and its cosmology. The equations of motion are described in section III. We consider cosmological 
era in section IV and numerical result are shown in section V. We conclude in section VI. 



II. EDDINGTON-BORN-INFELD COSMOLOGY 



In the EBI model studied here the action has three variables, the metric g plJ , the Born-Infeld connection C p pv and 
the matter field . The EBI action is 

%^,CV,*] - Y^G I dix {\J\^\ R + ^i2\/\9^- l2R ^\) + J d 4 z£m(*,. 9Al ,) (1) 

The action above has two extra constants, the length scale I which resembles the dimension of a length or 1/y/R, and 
a, which is a dimensionless parameter. The Born-Infeld Ricci tensor K^ v is symmetric under interchanging fj, and v 
which is a result of symmetric properties of the Born-Infeld connection C M „. As in standard general relativity, 

Kfn, — K p fipv (2) 

where 



kp — CP -X- CP C a — CP — CP C a d\ 

The conventional matter fields are included in the Lagrangian £ m . Since the two dynamical fields g^ v and C p ^ u are 
independent, the Born-Infeld connection can effectively be expressed in term of a new symmetric metric q lll/ (x), 

CV = \l pa (<^W + - - ( 4 ) 

giving a version of bi-metric theory. As in standard case, for the new metric covariant derivative vanishes 

D p q»v = 0, (5) 

where the covariant derivative is performed under the Born-Infeld connection, i.e. 

D p q^ v — d p q pv — C a pp q av — C° vp q ap . (6) 

Varying the action (1) with respect to two dynamical fields, the metric g pv and the connection C£„, yields the following 
equations of motion, 



Gv» = \l ^ S K( ^ 9„ ( -An?) " 9« v + S.GT™ . (7) 



Defining 



hence (7) can be written as 



G pv = -U^g^g^ + SnGT™. (9) 
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Varying of the action with respect to the connection C£ v , one can find D p (yfqq^ v ) = Taking determinant of (8) 
then we obtain, 

= ^(g^ + aq^) . (10) 

The first term in Eq.(9), is a modification from the Born-Infeld part. T™ is the matter field energy-momentum tensor. 
These results agree with the ones first shown in [9]. The two metrics g^ v and q pv possess homogeneity and isotropy 
with flat spatial curvature, 

g^dx^dx" = -dt 2 +a(t) 2 (dx 2 + dy 2 + dz 2 ), (11) 
q^dx^dx" = -X(t) 2 dt 2 +Y(t) 2 (dx 2 +dy 2 +dz 2 ) (12) 

g t t = —1, due to the gauge freedom in time. Here X(t) is the time rescaling of the metric q pu whereas a(t) and Y(t) 
behave like scale factors in g pv and q^ v respectively. The a(t ) is set to 1 so that H — a(t ) as in [9]. 

III. THE EQUATIONS OF MOTION 

Applying the metric ansatz Eqs. (11) and (12), to the equations of motion (9) and (10), we obtain first order 
equations, which are 



1 fY 3 \ 1 8ttG . 

H = )x + — < e ™ + !, - ) - < 13) 



dt\X J \Y 2 J 

X 2 ( 1 3a 2 \ 
[y) = W{ a -2X^ + 2Y^)- (15) 

It should be noted that I has dimensions of length (M whereas Y has dimensions of a. Finally, X is dimensionless. 
From Eq. (13), we can introduce an energy density as 

&E 8^G^' (16) 
and by taking the derivative of the Friedmann equation, one can find an effective pressure for this dark component as 

Px = w x gx, (17) 

where 

1 X' Y' 

and a prime denotes differentiation with respect to N = In a. 
Let us now introduce the variable 

_ 8nG gx 1 Y 3 
X ~ 3H 2 ^3l 2 H 2 Xa 3 ' [ ' 

In terms of this variable, the Friedmann equation can be written as 

1 = fix + + , (20) 

where we have defined, as usual, 

8TrGg m 8TrGg r 
m " 3H 2 ' = 3H 2 ' ^ ' 

and have assumed p m oc a~ 3 , p r oc a~ 4 . 
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We will demand Six > 0, as Six represents an effective matter density, otherwise Sl m ^ r could assume values larger 
than unity. From Eqs. (21), we find 



2H' 



0. 



Eq. (14) can then be rewritten as 



and Eq. (15) can be rewritten as 
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Therefore, we also need an equation for H. This can be found by differentiating the Friedmann equation as 

si'x + si' m + si' r = o, 



(22) 
(23) 

(24) 
(25) 



or 

2 jjt 

Sl' x = —{1- si x ) + 3(1 - Sl x ) + si r ■ 

Therefore the dynamical autonomous equations can be written as 

2H' 



SI'; 

SI' 



, - — (I -Six) +3(1 -Six) + Or, 
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(26) 

(27) 
(28) 

(29) 
(30) 



where we have introduced the dimensionless variable K 2 = Hfil 2 . This shows that the present value of H, can be 
re-absorbed into the free parameter K. In terms of these variables we find 



w x 



2 H' SI' 
-1 x 



3 H 3Sl x 



(31) 



IV. COSMOLOGICAL ERAS 



Let us consider the different eras in the cosmological history. We can distinguish the following cases. 
1. Radiation era: We can set Sl r = 1 and Six — 0. This fixes Sl m = 0. All the equations of motion are satisfied if 



H ~ 



which implies H = — 
2t 



(32) 



as expected. 
2. Matter era: Now we have two options: 
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(a) We can assume fl m — 1 and fl r = 0. In this case fix = 0. This implies that we are considering dark matter 
as an extra matter component (inside p m ) and not the X dark component. In this case the equations of 
motion are satisfied if 

^j=~\, that is ff=|, (33) 

as expected. 

(b) Now we assume that the dominant dark component behaves as dark matter, whereas fl m — > fib, that is, 
the matter component reduces to the baryon component and we suppose it is not the dominant one. In 
this case we need to impose fix = 1, and fl r = 0. Since we still want that H' /H = —3/2, the equations of 
motion cannot be solved at the same time. Therefore this case shows that if the X-components gives an 
effective dark-matter contribution in the past, it cannot be along a fixed point solution. However, there 
could be a transient solution from fix = and fix = 1 which could still mimic a dark-matter component. 

3. Dark energy era: In this case we set fl m = = fl r together with fix = 1. We look for de Sitter solution that is 
H' = 0. The equations of motion imply 

(34) 



H 2 K 2 



Therefore, X = X^s = constant. Then for a de Sitter solution we find 

-ds/- 

■A 



(35) 



Furthermore, the equations of motion give 



which implies that a < 1. We can also write 



H * s ~ 3K 2 (1 - a) ■ (37) 

4. Dark Energy for the case a > 1: Let us consider the case when, at very late times, X/H = A w constant, 
fix ~ 1, flr ~ w fl m , but still H 1 /H — > constant, as well as X'/X — > constant. Then, by neglecting any 
constant term with respect to the X term, we find the following two conditions which need to be satisfied 

3 1 / 3 2X' 



\i/3 K 4/3 X 

^ X' 



VTiXK X 
which imply 



— -3 = 0, (38) 
= 0, (39) 



3 ( 3 l/3 _ \*/3 K i/3) 2 

a = 4XV^ ■ (40) 



This solution is not a de Sitter solution, as in fact we find 



^-^- 32/3^/3^4/3 t- 1 - ( 41 ) 



V. NUMERICAL DISCUSSION 

Let us consider a numerical solution of the Eqs. (27)-(30). Rewriting Eq. (30) as 

and by allowing the constant / (or K) to take also negative values (but X^0), then we recover both the branches of 
Eq. (30). Notice that K < 0, Eq. (35), and implies X^s < 0, for a < 1. This further implies that the two branches, 
on their de Sitter solution, will differ by the sign of the final value of X. 
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A. Initial conditions 

Let us solve the equations of motion from a given redshift (z = Zi ~> 1), such that at z — Zi the universe is in the 
radiation era. We will set the initial condition for the Hubble parameter, during the radiation era, as the one given 
by GR, namely 

Hi = fff R) = yjsirfl + ^e-^ + (1 - fig*> - fi r , ) . (43) 

In what follows, we will fix the value of fl r ,i, during radiation era, at N = Ni = — log 10 (l + 1.76 x 10 5 ), such that 
Q r (N = 0) = £l r ,o which will be set equal to a fixed value. In this model, we have five parameters, f2 m ,o, Xi, £lx,i, 
a, K. However, we will fix the initial condition for Clx,i by requiring the condition Hk,o = 1 — ^m,o — ^r,o to hold. 
Finally, the four parameters, TOl o, Xi, a, K, will be considered to be free. In particular, since the A-component 
is supposed to explain both dark matter and dark energy, we will set the following range < O m .o < 0.4. We run 
the other parameters to change over a large range, —200 < Xi < 200, — 10 < a < 10, and — 15 < K < 15. Notice 
that in this parameter range, the system does not have a ACDM limit, therefore one expects deviations from the 
concordance model. Since this model has been introduced to explain dark energy and dark matter at the same time, 
this no-ACDM limit is in fact well motivated. 



B. Results 

We have calculated the total \ 2 f° r this model by using WMAP7 data (the background constraints on the two 
CMB shift parameters [23]), the BAO (SDSS7) data (two points) [24], and supernova type la (constitution data) [25], 
following the same method followed in [26]. The minimum for the x 2 is located at 

Q mfi = 0.250078 , K = 8.629636 , a = 2.760611 , X, = 77.73029 , where X 2 = Xmi„ = 484.505 , (44) 

where we have also fixed O ri j = 0.999827, and £lx,i — 1.14432 x 10 -6 for the reasons already explained above. Trying 
to set priors on f2 mi o like £l m ,o = Qb.o (i- e - fixing the scalar field to be the main source of dark matter) leads to much 
larger values for % 2 . Furthermore, data tends to prefer clearly the a > 1 case, as for < a < 1, the \ 2 increases. 

Nonetheless, the minimum value for xLin * s s * m much larger than ACDM's value (xacdm ~ 469). The x 2 f° r ACDM 
has two free degrees of freedom (fi TO o, A) that we can vary. Therefore according to the x 2 -probability distribution, at 
95% confidence level, ACDM rules out those models, at 2-cr, whose fit to the same data will lead to X 2 — Xacdm > 5-99. 
However, the models discussed here have \ 2 — 484.5, so that \ 2 ~ Xacdm = 15.5, which implies that these models 
are excluded at 2-a. This large difference implies the model under consideration does not fit the data, already at 2a, 
as well as ACDM. This is tantamount to saying that the ACDM cosmological evolution rules out this class of models. 
Since the \ 2 f° r the model studied here is higher than ACDM's one, we can deduce that data do not support well the 
evolution of the effective equation of state plotted in Fig. 2. It should also be pointed out that for the parameters 
for which x 2 = Xmim tne sca l ar field, although it has in the past wx ~ 0, it is anyhow a subdominant dark matter 
component (since on the minimum-x 2 solution, the dust-like dark matter contributes up to £l m ,o ~ 0.25). This implies 
that the scalar field starts dominating the evolution of the universe only at late times, that is it contributes to the 
dynamics essentially only as a dark energy field. But it is a dark energy field which, at early times, it is quite different 
from a cosmological constant: this may be part of the reason why, in this case, the model cannot not fit the data well. 
One option would be adding a bare cosmological constant (as also proposed in [16]), but in this case the model loses 
part of the interest as it would stop being an attractive dark energy model. Furthermore the evolution tends to lead 
to a fast transition of the effective equation-of-state parameter. This may also contribute to a worse fit to the data 
compared to ACDM. 

It should be noted that negative values for any of Xi, a, and K leads to very large values for x 2 (typically larger 
than 1000), giving a bad fit to the data. 

VI. CONCLUSIONS 

We have studied the Eddington-Born-Infeld scalar field which was proposed to model both dark matter and dark 
energy at the same time. We have solved the equations of motion and studied the behavior of the background at 
different times: at early times, indeed the scalar field behaves as a dark matter component with equation of state 
parameter wx ~ 0. Only at late times, the field can lead the dynamics of the universe to an accelerated regime, which 
depending on the parameters of the model, is described by either a de Sitter solution, or a rather different dynamics 
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Figure 1: Evolution for H'/H = H/H 2 (left panel), and for the variable X' /X (right panel). The evolution, starting from 
radiation domination, passing through matter domination, at late times, tends to a super-accelerating final state on the 
minimum-x 2 solution. 



W X 




N N 

Figure 2: Effective equation of state for the scalar field (left panel), and plot of the ratio X/H (right panel). At early times the 
Eddington-Born-Infeld scalar field behaves as a dark matter component (wx ~ 0), then, at late times, it drives the evolution 
of the universe. This plot shows the evolution for the parameters which minimize the \ 2 given in (44). Notice that since 
a > 1, the final state is not de Sitter (wx < —1), rather it tends to the solution characterized by H'/H — > constant, and 
X/H — > constant (right panel). 



described by H/X — > constant, and H'/H — > constant (where a prime denotes differentiation with respect to the 
JV = Ina). 

The fact that, at early times, the scalar field behaves as a dust component can in principle alleviate the problem 
of finding a dark matter component, as indeed the nature of dark matter and dark energy would have the same 
explanation. 

In order to see whether this model is viable or not, we studied the cosmological constraints that its dynamics have 
to pass when considering WMAP7 data, BAO and Supernova type la. For this goal, we have calculated the x 2 as 
a function of four free parameters, that is K 7 a (two theoretical dimensionless parameters of the model) together 
with Cl m ,o (which states how much of an extra standard dust component is needed), and Xi, the initial value for the 
time-rescaling component of the Born-Infeld metric. 

We have found that the model cannot give a good fit to the data (compared to ACDM), and hence the model 
cannot be considered viable. We have proved this statement by constraining the background. This approach differs 
from the one followed in [16], where the authors studied the evolution of the cosmological perturbations in order to 
constrain the growth of structures. In particular, we have used only the constraints on the background coming from 



8 



the WMAP7 data. Instead in [16], the authors considered constraints only on the perturbations power-spectrum. It 
is possible, as also suggested in [16] that introducing a cosmological constant would improve the fit, but, on the other 
hand, the model would partially lose its original motivation of explaining at the same time both dark energy and dark 
matter. In particular data prefers the non-de Sitter solution, preferring a fast transition to values for wx < — 1- In 
this model, at early times, a cosmological constant is absent from the beginning, as the scalar field initially (and up 
to very recently) behaves as a dark matter component. 
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